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Features of two-pion Bose-Einstein correlations based on event-by-event analysis in 

smoothed particle hydrodynamics 
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We examine the space-time structure of the particle-emitting sources with fluctuating initial 
conditions in smoothed particle hydrodynamics. The two-pion correlation functions of single events 
for the sources exhibit event-by-event fluctuations. The large event-by-event fluctuations and wide 
distributions of the error-inverse-weighted fluctuations between the HBT correlation functions of 
single and mixed events are important features for the sources with the fluctuating initial conditions. 
The root-mean-square of the weighted fluctuations is a signal to detect the inhomogeneity of the 
systems produced in high energy heavy ion collisions. 
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I. INTRODUCTION 

The main physics goal of heavy ion collisions at the 
BNL Relativistic Heavy Ion Collider (RHIC) is the study 
of the extremely hot and dense matter, the quark-gluon 
plasma (QGP), formed in the early stage of the colli- 
sions. Hydrodynamics may provide a direct link between 
the early state of the matter and final observables and 
has been extensively used in high energy heavy ion colli- 
sions. The hydrodynamic calculations with the equation 
of state of the QGP agree well with the RHIC t>2 data 
of the elliptic flow at low transverse momentum p T < 2 
GeV [l|, S S 0, Q , which is believed as an important ev- 
idence of the existence of a strongly coupled QGP in the 
early stages of the collisions [sHdii&E]- However, hy- 
drodynamic results can not explain the RHIC Hanbury- 
Brown-Twiss (HBT) measurements, Rout/Rside ~ 1 
[HE EH- R is the so-called HBT puzzle. 

In Ref. E a granular source model of QGP droplets 
evolving hydrodynamically was put forth to explain the 
HBT puzzle. The suggestion was based on the obser- 
vation that in the hydrodynamic calculations for the 
granular source the average particle emission time scales 
with the initial radius of the droplet, whereas the spa- 
tial size of the source is the scale of the distribution of 
the droplets. For a granular source with many of the 
small droplets distributed in a relatively large region, 
the HBT radius i? ou t can be close to i? s ide El- In Ref. 
[lq the authors further investigated the elliptic flow and 
HBT radii as a function of the particle transverse mo- 
mentum for an improved granular source model of QGP 
droplets. They argued that although a granular struc- 
ture was suggested earlier as the signature of a first-order 



phase transition E E E E E E E E HI, the 
occurrence of granular structure may not be limited to 
first-order phase transition [16, 25]. In an event- by- event 
basis, the initial transverse distribution of energy den- 
sity in nucleus-nucleus collisions have been known to be 
highly fluctuating EES1- This lar S e spatial fluctu- 
ation may facilitate the occurrence of instability of the 
system during it violent expansion subsequently and the 
fragmentation to ma ny gran ular droplets together with 
surface tension effect [16!, 25]. So, the examination of the 
space-time structure in event- by- event basis is important 
to understand the system evolution and the HBT puzzle. 

Smoothed Particle Hydrodynamics (SPH) is a suitable 
candidate that can be used to treat the system evolu- 
tion with large fluctuating initial conditions for inves- 
tigating event-by-event attributes in high energy heavy 
ion collisions [23, [29j]. The main idea of SPH is the 
parametrization of the fluid in terms of discrete La- 
grangian coordinates attached to small volumes (call 
'particles') with some conserved quantities in hydrody- 
namics m, [H, El. NEXSPHERIO is a SPH code 
[H, SO, 



28, [291, 130, 

321 ] with event-by-event initial conditions gen- 
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erated by NEXUS event simulator [27|. It has been used 
to study a wide range of problems in high energy heavy 
ion collisions E SS HE tH E H El ■ In this article 
we use NEXSPHERIO to simulate the evolution of the 
system produced in the collisions of ^/snn = 200 GeV 
Au+Au at RHIC. We will examine the space-time struc- 
ture of the systems and investigate the two-pion HBT 
correlation functions in event-by-event basis. Our re- 
sults show that the systems are inhomogeneous in space 
and time both for a first-order and a cross-over transi- 
tions between the QGP and hadronic phase. There are 
"lumps" in different regions of the system. The two-pion 
HBT correlation functions in event-by-event basis exhibit 
large fluctuations, which lead to a wide distribution of 
the error-inverse- weighted fluctuations between the HBT 



correlation functions of single and mixed events. These 
large fluctuations and wide distributions are important 
features for the sources with event-by-event fluctuating 
initial conditions. The root-mean-square of the weighted 
fluctuations is a signal to detect the inhomogeneity of the 
systems produced in high energy heavy ion collisions. 



II. SYSTEM EVOLUTION AND SPACE-TIME 
STRUCTURE 

In hydrodynamics the behavior of system evolution is 
determined by the initial conditions and the equation of 
state (EOS) of the system. The system initial states of 
NEXSPHERIO are given by the NEXUS code, which can 
provide detailed space distributions of energy-momentum 
tensor, baryon-number, strangeness, and charge densi- 
ties at a given initial time (which is taken to be Tn = 1 
fm in our calculations), in event-by-event basis [271. I28J]. 
In our calculations we use two kinds of EOS. The first 
one, EOS-I, considers a first-order phase transition at 
T c = 160 MeV between a QGP and a hadronic resonance 
gas as used in Refs. [33, |33|. The QGP is an ideal gas 
of massless quarks (u, d, s) and gluons and the hadronic 
gas contains the resonances with mass below 2.5 GeV/c 2 , 
where volume correction is taken into account [3ll I33J ] . 
Because at RIIIC energy the net baryon constituent is 
very small and the phase transition is a smooth cross 
over, in the mid-rapidity region of the heavy ion colli- 
sions, a modification for the EOS with first-order phase 
transition in the hydrodynamical code should be consid- 
ered |28l . [34l | . Accordingly, we introduce EOS-II which 
is obtained by smoothing the EOS-I in the transition re- 
gion with the entropy density suggested by QCD lattice 
results [15|, [3y, [33, [38|, [39| ■ Fig. 1 shows the relations of 
the thermodynamical quantities for the systems evolving 
with the EOS-I (grey) and the EOS-II (black), where s, e, 
and P are entropy density, energy density, and pressure 
of the system. The width of the transition AT for the 
EOS-II is taken as 0.1T C [1|. 

The coordinates used in NEXSPHERIO are r = 
z 2 , x, y, and r) = (1/2) ln[(t + z)/{t - z)} [H, [H, 



VT 2 r 

\3l\ . They are convenient for a system with rapid lon- 
gitudinal expansion. However, in order to examine the 
space-time structure of the system with nonlocal coordi- 
nates we work in the center-of-mass frame of the system. 
Figs. 2(a), (b), (c), and (d) are the pictures of the trans- 
verse distributions of energy density for one NEXSPHE- 
RIO event (^/sjvjv = 200 GeV Au+Au) with impact pa- 
rameter b — f m and evolving with EOS-I. The pictures 
are taken for the region (\x,y\ < 12 fm and \z\ < 1 fm) 
and at t = 1, 5, 9, and 13 fm/c with an exposure of 
At = 0.3 fm/c. Figs. 2(a'), (b'), (c'), and (d') are the 
pictures for the event with the same initial conditions 
but evolving with EOS-II. One can see that the systems 
are inhomogeneous in space and time both for the events 
evolving with the EOS-I and the EOS-II. There are many 
"lumps" in the systems. By comparing the pictures for 
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FIG. 1: (a) entropy density s, (b) energy density e, (c) pres- 
sure P, and (d) P/t of the systems evolving with EOS-I (grey) 
and EOS-II (black). 



different events we find that the lump locations arc differ- 
ent event-by-event. Figs. 3((a)-(d)) and ((a')-(d')) show 
the pictures of the transverse distributions of energy den- 
sity for the events with different impact parameters and 
evolving with EOS-I and EOS-II, respectively. One can 
see that the number of the lumps in the system decreases 
with impact parameter increasing. Because the systems 
evolving with the two kinds of EOSs have the similar 
space-time structure, we will consider only the EOS-II in 
our calculations later. 



III. TWO-PION CORRELATION FUNCTIONS 
IN EVENT-BY-EVENT BASIS 

Two-pion interferometry is a powerful tool for prob- 
ing the space-time structure of particle-emitting source. 
For the source with dense and void density oscillations, 
the single-event HBT correlation functions will appear 
fluctuations relative to the mixed-event HBT correlation 
functions [23, [24[ . The fluctuation patterns are different 
event- by-event. 

Assuming that final identical pions are emitted at a 
system configuration characterized by a freeze-out tem- 
perature Xf, we may generate the pion momenta ac- 
cording to Bose-Einstein distribution and construct the 
single-event and mixed-event two-pion correlation func- 
tions [H, Si] for NEXSPHERIO events. Fig. 4 shows 
the two-pion correlation functions C(q S idc, <Zout,9iong) for 
the NEXSPHERIO events with different impact param- 
eters. Here <7 s idc, <2out, and QWg are the components of 
"side" , "out" , and "long" of relative momentum of pion 
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FIG. 2: Transverse distributions of energy density of the 
events evolving with EOS-I ((a)-(d)) and EOS-II ((a')-(d')) 
at different times for ^/snn = 200 GeV Au+Au collisions, 
6 = fm. The unit of energy density is GeV/fm . 



FIG. 3: Transverse distributions of energy density of the 
events for ^/sjvjv = 200 GeV Au+Au collisions with 6=5 
and 10 fm. The systems evolve with EOS-I ((a)-(d)) and 
EOS-II ((a')-(d'))- The unit of energy density is GeV/fm 3 . 



pai r [401 . |41j . which are calculated in the LCMS frame 
[13 |. In each panel of Fig. 4, the dashed lines give 
the correlation functions for a sample of different sin- 
gle events and the solid line is for the mixed event ob- 
tained by averaging over 40 single events. In the present 
calculations, the freeze-out temperature is taken to be 
150 MeV. For each single event the total number of 
generated pion pairs in the relative momentum region 
(q s idc,qout,qiong < 200 MeV/c) is N^ = 10 6 and the 
numbers of the pion pairs in the relative momentum re- 
gions (qt < 200MeV/c; q v q k < 30MeV/c) are about 
2.7%Nn Tn where i, j, and k denote "side", "out", and 
"long". 

From Fig. 4 it can be seen that the correlation func- 
tions for single events exhibit fluctuations relative to 
those for mixed events. These fluctuations are larger for 
bigger impact parameter. It is because that the number 
of the lumps in system decreases with the impact param- 
eter and the fluctuations are larger for the source with 
smaller number of lumps [23l . |24|. It also can be seen 
that in the longitudinal direction the correlation func- 
tions exhibit oscillations which can not be smoothed out 
by event mixing. This is because that there are two sub- 
sources moving forward and backward the beam direc- 
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FIG. 4: Two-pion correlation functions for a sample of differ- 
ent single events (dashed lines) and mixed events (solid lines) 
for the NEXSPHERIO events with different impact parame- 
ters. 
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FIG. 5: Two-pion correlation functions for a sample of differ- 
ent single events (dashed lines) and mixed events (solid lines) 
with the initial rapidity cut r/o > 0. 



tion. By applying an additional cut for the initial rapid- 
ity of the "smoothed particles" , r/ > 0, we find that the 
oscillations of the mixed-event correlation functions are 
smoothed out as shown in Fig. 5(c). 



IV. DISTRIBUTION OF THE FLUCTUATION 

OF SINGLE-EVENT HBT CORRELATION 

FUNCTION 

In last section we show that the two-pion correlation 
functions of single events exhibit event-by-event fluctua- 
tions. In realistic experimental situations, the number of 
identical pion pairs in each single event is limited. Tra- 
ditional HBT measurements are based on mixed-event 
analysis because of statistics. The event-by-event fluc- 
tuations are smoothed out in the mixed-event analysis. 
In order to observe the event-by-event fluctuations, we 
next investigate the distribution dN/df of the fluctu- 
ations between the correlation functions of single and 
mixed events. | C s (<?j) — C m (qi)\, with their error- inverses 
as weights [24j . 



/(ft) 



\C s (qj) - C m (qi)\ 
A|C a ((fc)-C m (ft)| 



(1) 



In the calculations, we take the width of the relative mo- 
mentum qi bin as 10 MeV and use the bins in the region 
20 < qi < 200 MeV. The up panels of Fig. 6 show the 
distributions of / in the "side" , "out" , and "long" di- 
rections, obtained from 40 simulated vents. The impact 
parameter for these events is 5 fm (b = 5 fm) and the 
number of correlated pion pairs for each of these events 
is 10 7 (Nmr — 10 7 ). The thick lines are the results for the 
events with fluctuating initial conditions (FIC). Because 
in the last analysis the fluctuations of the single-event 
correlation functions are from the FIC, for comparison 
we also investigate the distribution dN/df for the events 
with the smooth initial conditions (SIC) obtained by av- 
eraging over 30 random NEXUS events [28|, [33, [34| (al- 
though it is not a realistic case) . It can be seen that the 
/ distributions for FIC are wider for the higher freeze- 
out temperature T{ — 150 MeV than those for the lower 
T{ = 125 MeV, and both for the freeze-out temperatures 
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FIG. 6: The distributions dN/f for 40 events with FIC and 
SIC. 6 = 5 fm. 
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FIG. 7: The distributions dN/df for 40 events, where / are 
calculated for the variables of (/trans and q. b = 5 fm. 



the distributions for FIC are much wider than the cor- 
responding results for SIC. The down panels of Fig. 6 
show the distributions of / for 40 simulated events with 
b = 5 fm and N nir = 5x 10 6 . It can bee seen that the 
distributions for FIC are also wider than those for SIC 
in this case. 

In experiments the number of correlated pion pairs in 
one event, N vn , is limited. It is related to the energy 
y/SNN of the collisions. For a finite iV T7r , sometimes 
we have to reduce variable numbers in data analyses al- 
though it will lose some details. In Fig. 7 we show the dis- 
tributions of / for the variables of the transverse relative 
momentum gtrans and the relative momentum q for 40 
events with b = 5 fm. One can see that for N n „ = 5 x 10 6 , 
the distributions for FIC are much wider than those for 
SIC both for qtrans and q. Even for N„ n = 5 x 10 5 , the 
widths for FIC arc visibly larger than those for SIC for q. 
In order to examine quantitatively, we calculate the root- 
mean-square (RMS) of the /. Figure 8 shows the RMS 



E4 
2 


2 
il-5 



^ 



1 
0.5 



(a) for q 
- b=9 fm 




A -A A 

~""l I 



■(c) for q 
' b=5 fm 




i i i i i "in 



( b ) for <7tra„s 
- b=9 fm 

G— © FIC 
A- ^ SIC 




--A-A A 

"mil I I'm 



(d) for q t 
v ' v trans 

" b=5 fm 




' ' i i ' "mi 



1 10 50 1 10 50 

A^(xl0 5 ) A^(xl0 5 ) 



FIG. 8: The root-mean-square of / for the variables of q and 
<?trans as a function of N-xn- Tf = 150 MeV. 



/ rms as a function of N vv for the 40 events with 6 = 9 
fm (the up panels) and 6 = 5 f m (the down panels). The 
freeze-out temperature is T{ — 150 MeV and the SIC are 
obtained by averaging over 100 NEXUS events. It can 
been seen that the values of / rms rapidly increase with 
N^TT for FIC because the errors in Eq. (JTJ) decrease with 
N V7r . For FIC the results for b — 9 fm are larger than 
those for b = 5 fm correspondingly because the differ- 
ences \C s {qi) — C m (qi)\ in Eq. (Q]) increase with b. For 
SIC the values of / rms are almost independent from N„ n 
. It is because that both the differences and their errors 
in Eq. |T]) decrease with N W7r in the SIC case. 

At RHIC energy the event multiplicity of identical pi- 
ons, M w , is about several hundreds for central collisions. 
The order of N^ is about 10 5 (~ Ml/2). However, at 
the higher energy of the Large Hadron Collider (LHC) 



at CERN, M v will be about two thousands and the or- 
der of -/V^tt will be 10 6 . Our results indicate that if the 
particle-emitting sources at LHC energy have the similar 
inhomogeneity as that at RHIC energy, the RMS of / at 
LHC energy will be much larger than the corresponding 
result at RHIC energy. The distributions of / and their 
RMS are useful observables for the inhomogeneity of the 



V. CONCLUSION 

In an event-by-event basis, the initial density distri- 
bution of matter in nucleus- nucleus collisions have been 
known to be highly fluctuating. The fluctuating ini- 
tial conditions lead to event-by-event inhomogeneous 
particle-emitting sources. For these sources the two-pion 
HBT correlation functions of single events exhibit event- 
by-event fluctuations. Because of data statistics tradi- 
tional HBT measurements are mixed-event analysis. In 
these measurements the fluctuations of the correlation 
functions of single events are smoothed out. In order to 
observe the fluctuations we examine the distributions of 
the error-inverse- weighted fluctuations between the HBT 
correlation functions of single and mixed events. We find 
that the distributions and the RMS of the weighted fluc- 
tuations are useful observables for the inhomogeneity of 
the sources. If the particle-emitting sources produced 
in the heavy ion collisions at LHC energy are inhomo- 
geneous, these RMS will be much larger than those at 
RHIC energy. 
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